Orbital coupling and superconductivity in the iron pnictides 
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We demonstrate that strong inter-orbital interaction is very efficient to achieve superconductivity 
due to magnetic fluctuations in the iron pnictides. Fermi surface states that are coupled by the 
antiferromagnetic wave vector are often of different orbital nature, causing pair-hopping interactions 
between distinct Fe-3d orbitals to become important. Performing a self-consistent FLEX calculation 
below Tc we determine the superconducting order parameter as function of intra- and inter-orbital 
couplings. We find an -pairing state with Tc — 80K for realistic parameters. 



High superconducting transition temperatures and the 
proximity to antiferromagnetic order strongly suggest an 
electronic pairing mechanism in the FeAs systems|[l|. The 
vicinity to a spin density wave instability with param- 
agnons as dominant collective mode is key for spin fluc- 
tuation induced superconductivity, where pairing is the 
result of paramagnon exchange. However, superconduc- 
tivity in the iron pnictides occurs not only in the im- 
mediate vicinity of the magnetically ordered state and 
the viability of spin fluctuation induced pairing becomes 
an issue that requires a quantitative analysis. In ad- 
dition, multi-orbital effects of the Fe-Sd bands with a 
fining of approximately six electrons per Fe-site add to 
the complexity of these systems: Electronic structure 
calculations 0, 0| yield two sets of Fermi surface sheets, 
one around the center of the Brillouin zone (F-point) and 
the other around the M-point, shifted from F by the 
magnetic ordering vector Qj3l- Inter-band scattering of 
electrons has been proposed to lead to unconventional 
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parameters other solutions exist]; 
coupling tends to support the -pairing state where the 
gap functions on the two Fermi surface sheets have op- 
posite sign. 

Crucial for all scenarios based upon inter-band scatter- 
ing is that states IV'r.k) on one Fermi surface are coupled 
to states IV'M.k+Q) on another Fermi surface, and vice 
versa. The natural starting point to describe electron- 
electron interactions in transition metals is however not 
in terms of bands, but rather in terms of local orbitals 
\a). Here a — xz, yz, xy, x^ — y^ and — refers 
to the Fe-Sd orbitals, with intra- and inter-orbital direct 
Coulomb interactions, U and U' , as well as Hund's rule 
coupHng Jh and inter-orbital pair hopping J'. The im- 
portance of orbital effects in the iron pnictides was also 
stressed in Ref . (l3| . As we will see below, the dominant 
effective spin-fluctuation induced pairing interaction in a 
multi-orbital system is of the pair-hopping form: 
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ing interactions, while a b corresponds to an inter- 
orbital pairing interaction. In both cases. Cooper pairs 
are predominantly made up of electrons in the same or- 
bital: {d-kaidkai) 7^ 0. In the band picture this yields 
the inter-band pairing interaction 
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A pair of electrons in orbital b is scattered into a pair 
in orbital a. For a = b we consider intra-orbital pair- 



The dominant momentum transfer in the spin fluctua- 
tion approach is of course k — k' = Q. It is interesting to 
observe that electronic structure calculations show that 
IV'r.k) and IV'M.k+q) are often dominated by different 
orbitals. For example, if {x z\^|JT ,'k) is large, it holds that 
{xz\ipM,'k+Q) for the same k is small, while {xy\ipAi,-k+Q) 
might be large. In Fig. [l^ we illustrate this effect where 
distinct colors refer to the dominant orbitals on the Fermi 
surface. We used the tight binding parametrization of 
the flve band model of Ref.[l3|, where a similar plot was 
presented. The three dominant orbitals on the Fermi sur- 
face are xz, yz, and xy. Connecting a Fermi surface point 
by Q = (tJ", 0) or (0, tt) leads in most cases to a different 
orbital. Thus, the orbital composition of the wave func- 
tion at the Fermi surface frustrates intra-orbital pairing. 
In other words, inter-band but intra-orbital scattering of 
spin fluctuations (oc W^^., ) provides a less efficient pair- 
ing glue, if compared to inter-orbital scattering cx W^\_, 
{a ^ b) of equal size. It is crucial to determine un- 
der what conditions collective paramagnons with strong 
inter-orbital pair-hopping exist. 

In this Rapid Communication we solve the two-orbital 
many-body problem in the superconducting state for 
varying intra-orbital (U) and inter-orbital {U' , J = Jh = 
J') couplings within the self-consistent fluctuation ex- 
change (FLEX) approximation (lil . We obtain supercon- 
ductivity with s*-pairing. The superconducting order 
parameter is determined self-consistently and vanishes at 
Tc ~ 80K. We demonstrate that strong collective inter- 
orbital spin-fluctuations are efficient to increase super- 
conductivity. To solve the FLEX equation on the imagi- 
nary frequency axis for a lattice of iV = 32 x 32 sites and 
at temperatures as low as T ~ lOK we require 2^^ = 8192 
Matsubara frequencies. At the moment, this restricts our 
analysis to consider only two orbitals. In Fig. [TJj we show 
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Figure 1: (Color online) Fermi surface of the five-orbital (a) 
and two-orbital (b) tight binding model of the Fe-3ci states 
in the unfolded Brillouin zone (one iron atom per unit cell). 
Colors indicate the dominant orbital that contributes to the 
bands: xz (orange/gray), yz (blue/dark gray), and xy (yel- 
low/light gray). The antiferromagnetic vector Q — (7r,0) 
mostly connects states dominated by different orbitals. The 
tight binding parameters are from Ref.[13] for panel (a) and 
Ref.[i3| for panel (b). 



the Fermi surface of a two band model with d^z and d, 
orbitals. The mentioned frustration of intra-orbital pair- 
ing is less pronounced for this simplified model. Yet, the 
phase space for inter-orbital pairing interactions is still 
larger compared to intra-orbital interactions. 
The model: We consider the Hamiltonian 

H = ^ Skdl^^di^ba - Jh ^ (^Sia ■ Sib + ^nian.ibj 
k.a6,cr i.a>b 

+ U ^^Uia^riiai + U' riiaiMb + J' bj^bib, (3) 
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where riiaa = dl^^diaa is the occupation of the or- 
bital a with spin a at site i and riia = '^^niaa- 
Sta = \ Yl,aa' diaa'^cra' diaa' is the electron Spin and bl^ = 
'^lai^iai the pair creation operator, respectively. For 
the tight binding band structure we use — — 
— 4:t4 sin kx sin ky and e^^ = — 2tiCosfca — 2t2COsfca- — 
4^3 COS kx cos ky — /i, where a — x (y) stands for xz (yz) 
orbital as well as the momentum coordinate with x — y 
and y — X. We use ti = —0.33 eV, = 0.385 eV, 
ts = -0.234 eV, and = -0.26 eV of Ref. Our 
results were obtained for a filling of n = 1.88 electrons 
per site, corresponding to moderate hole doping. The 
filling of a subset of bands is primarily determined to 
reproduce realistic Fermi surface geometries and yields 
commensurate magnetic fiuctuations. 

The multi- orbital fluctuation exchange approach: The 
FLEX equations for a multi-orbital problem are given in 
Ref. [l^ ■ In the normal state one obtains the single par- 
ticle self energy S^^ which yields the single particle prop- 
agator G^^. Here k = (k,zti;„) stands for the momentum 
vector k and the Matsubara frequency ujn = (2n -|- 1) ttT 
with temperature T. As it is important for our subse- 



quent discussion we summarize the key equations that 
occur in the superconducting state and determine the 
anomalous self energy: 
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This equation is the strong coupling version of the gap- 
equation. . . . — ^. . . stands for the summation 
over momenta and Matsubara frequencies. F^^ is the 
anomalous Green's function, that determines the Cooper 
pair expectation value (dkaTt^-khi) = Tj^n^k^- 
thermore, Y'^'^''^'' is the dynamic pairing interaction that 
depends on momentum, frequency and the orbital states 
involved, where q = (q,w„) with Vn = 2n7rT. Intro- 
ducing the two-particle quantum numbers A — (a, c) and 
B = {d, b) that label the rows and columns of two-particle 
states, the interaction, F^'^ = pac,db ^^ggQjjjes a m? x m?- 

dimensional symmetric operator F^, where m is the num- 
ber of orbitals. It is now straightforward to sum particle 
hole ladder and bubble diagrams. It follows: 
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where interactions in the spin and charge channel are: 

^s{c),q — '^Us(c) (y^T Xs(c),qUs(c)j Xs{c),qUs(c) 

-\Us{c){Xs,q-Xc,q)U,(^c)- (6) 

Us and Uc are also m? x w'^-dimensional matrices of 
the interaction in the spin and charge channel, respec- 
tively. Close to a magnetic instability, the dominant 
contribution to F^ comes from the spin channel due 



to the Stoner enhancement ( 1 — Xs,qUi 



(l-Xs,qUs) 



For the in- 



teraction matrix in the spin sector holds [/°°^°° = JJ , 
while {or a ^ b holds Uf'''^ = U' , Uf^'"' = J' and 
IjaaM ^ rjj^g Hartree-Fock contribution Vhf ^ 

{Us + Uc)/2 is suppressing superconductivity, an effect 
caused by the repulsive direct Coulomb interaction. We 
find that the impact of direct Coulomb interaction is 
strongly reduced in the s"^ state with small, but finite, 
average {dial dial) fo'^ local Cooper pairing. For a dis- 
cussion of this Coulomb avoidance see Ref. Fi- 
nally, the irreducible particle hole bubble Xs(c).q is de- 
termined by normal and anomalous Green's functions: 
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reversal invariance and singlet pairing. We solved the 
set of coupled FLEX equations self consistently in the 
superconducting state. 

The pairing state and its T-dependence: The momen- 
tum dependence of the anomalous self energy is 
shown in the insets of Figs. [3] and [H The symmetry 
of <i>^^ and of the Hamiltonian are the same, correspond- 
ing to s-wave pairing[l3|. Nevertheless, the sign of is 
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Figure 2: (Color online) Temperature dependence of the 
anomalous self energy $£=o,7rT: proportional to the super- 
conducting order parameter, for U = 1.5eV, U' = 1.2eV, and 
J = 0.8eV . The inset shows the increase of with increas- 
ing inter-orbital coupling J at T = 60K. 



opposite on Fermi surface sheets around T and M, i.e., we 
obtain the pairing state that was proposed in Ref. 01 ■ 
While in general such a state can have nodes of the gap 
on the Fermi surface, our solution corresponds to a fully 
gapped state. For a recent discussion of the s^-state see 
Ref. In Fig. [2] we show the temperature dependence 
of the anomalous self energy ^IZo^ which is proportional 
to the superconducting order parameter. The feedback 
of the opening of a pairing gap onto the dynamic pairing 
interaction leads to the rather rapid growth of the order 
parameter below Tcfl^l- Tc ~ 80K is indeed of the correct 
order of magnitude. In the normal state the dynamics 
of paramagnons is overdamped Tq'^"" (1 + |a;|/a;s) ^ . 
For the parameters of Fig. [2] we find Ug {Tc) — 37meV. 
This energy scale is reduced compared to the typical elec- 
tronic energies because of the Stoner enhancement. It 
sets the scale for the Lorentzian lineshape of inelastic 
neutron scattering at Q above Tc. A strong pairing in- 
teraction is caused by a significant Stoner enhancement, 
which is controlled by the value of the magnetic inter- 
orbital coupling, J. This is demonstrated in the inset of 
Fig. [2] where we show the sensitivity of superconducting 
order with respect to J at T = 60K. 

Intra- versus inter-orbital pairing: For the pnictides, 
the pairing vertex Tq in Eq. ^ is dominated by 
a few matrix elements. We find that xf,',^ at q ~ Q 
has comparable diagonal elements Xd and somewhat 
smaller counter diagonal elements Xd two-particle 
space, while all other matrix elements are negligible. 
It then follows from Eq. ^ that the dominant ma- 



trix elements of Tq are F^'''''". 



If one interprets F^^'^" 
as effective low energy interaction, the combination of 
orbital indices yields precisely the pair-hopping form 
Eq. (J]) with W^i, = Ft'^t?. The effective Stoner 




Figure 3: (Color online) Intra-orbital and inter-orbital pair- 
ing interaction at zero frequency as function of momentum 
for U = 1.5eV, U' = leV and J = leV at a temperature of 
T = 70K along with the self-consistently determined anoma- 
lous self energy $£^.^t (inset) that determines the sign and 
momentum dependence of the superconducting gap. The 
pairing interaction is peaked for momenta Q = (tt, 0) and 
(0, tt). The anomalous self energy corresponds to -pairing 
with opposite sign of the gap on Fermi surface sheets around 
F and M. Inter-orbital pairing is significantly smaller. 
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enhancements for same orbitals a 
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{U + JH){l-{U + JH){xd + Xi)) 
follows Tf'- ^ iU' + J') (1 - ([/' + J') {xd + Xd)) 
Whenever U is significantly larger than the other cou- 
plings, Eq. im is dominated by interactions within 
the same orbital. As mentioned, the intra-orbital pair- 
ing interaction is however rather inefficient. The sit- 
uation changes when we consider comparable values 
for the intra- and inter-orbital Coulomb enhancements: 
U + Jh ^ U' + J', i.e. a regime with strong orbital fiuctu- 
ations. The pairing interaction is enhanced, as the nature 
of the wave functions on the Fermi surface can efficiently 
take advantage of coupling between distinct orbitals, see 
Eq. ^ and Fig.l. 

The condition U + Jh ^ U' + J' is at variance with 
the relations U — U' + 2J and J' — Jr that re- 
sult from the rotational symmetry of the bare Coulomb 
interaction[l^, if combined with evidence for sizable 
Hund coupling 13, 21 1 ■ We stress however that the inter- 
action parameters that enter an approximate theory such 
as FLEX are not identical to the bare Coulomb matrix 
elements [i^l- FLEX ignores crucial vertex corrections 
and U , U' and J should rather be considered low energy 
interaction parameters that have been renormalized by 
high energy excitations. Performing a renormalization of 
the Coulomb interactions within a multiband version of 
the Kanamori scattering matrix approach[23j, we indeed 
find that U + Jh — U' + J' for realistic values of the bare 
Coulomb matrix elements of FefiSl- In this approach, 
particle-particle excitations couple states with k and — k, 
i.e., states of same orbital nature. This reduces U more 
strongly than U', Jh, and J' . Thus, constraints due to 
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Figure 4: (Color online) Same as Fig. [3] but for parameters 
U = U' = 1.5eV and J = leV with same temperature T = 
70K. Now intra- and inter-orbital couplings are of comparable 
size, leading to the same superconducting order parameter 
as for Fig. [S] however for much smaller pairing strength, 
associated with moderate antiferromagnetic fluctuations (see 
the different scale compared to Fig. [3j. 

rotational invariance do not apply for low energy vertices 
that enter FLEX, and intra-orbital and inter-orbital pair- 
ing interactions can easily be comparable. The underly- 
ing low density expansion makes the Kanamori scattering 
matrix renormalization a very sensible approach for the 
pnictides, given their near semimetallic electronic struc- 
ture. 

In Figs. [3] and [4] we compare two different parameter 
sets that yield almost the same value for ^JJ" at the same 
temperature. The first case has a predominant intra- 
orbital interaction which needs to be very large in order 
to achieve pairing. In the other case the same pairing 
amplitude is obtained from intra- and inter-orbital pair- 
ing interactions. r°^''"' in Fig. |4] are almost one fifth of 
the pure intra-band interaction in Fig. [3l demonstrating 
the efficient role played by inter-orbital magnetic pairing 
interactions in the iron pnictides. 

In summary we presented a self-consistent FLEX anal- 
ysis of a two orbital model of the FeAs systems in the 
superconducting state. We determined the temperature 
dependence of the superconducting order parameter and 
showed that Tc — 60 — 80K, of the order of the experi- 
mental values, are clearly possible. The pairing state is 
with opposite sign of the gap on Fermi surface sheets 
around F and M|3|. In the iron pnictides, states that 
are coupled by the antiferromagnetic wave vector are 
often dominated by different local Fe-3d orbitals. This 
makes a purely intra-orbital pairing interaction quite in- 
efficient. Inter-orbital pairing due to antiferromagnetic 
fiuctuations yields the same pairing amplitude for much 
smaller Stoner enhancement, i.e., for more moderate val- 
ues of the magnetic correlation length. We expect this 
effect to be even stronger in a more realistic five-orbital 
description of the iron pnictides. Collective low energy 



pairing interaction between like and unlike orbitals, i.e., 
strong orbital fiuctuations, significantly enhances the vi- 
ability of the spin fiuctuation approach for superconduc- 
tivity in the pnictides. 
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